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ABSTRACT 

This  paper  focuses  on  the  effects  of  chemical  reaction  and  radiation  absorption,  on  unsteady  MHD  free 
convection  fluid  flow,  embedded  in  a  porous  medium  with  time  dependent  suction  and  temperature  gradient  heat  source, 
in  two  different  cases  of  boundary  conditions. 

Case  (I):  U  -  U  p  ,  0=\  +  Ee'",  C  =  l  +  &>'"  at  y  =  0 


u  — ^  0 , 

0^0, 
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as  y  — > 

Case  (II):  U  =  1  +  £B 

0  = 

\  +  eent , 

C  =  \  +  sen' 

u  — ^  0 , 

0^0, 

C^O 

as  y  — >< 

The  parabolic  partial  differential  equations  governing  the  fluid  flow,  heat  transfer  and  mass  transfer  have  been 
solved,  by  using  the  perturbation  technique.  The  velocity,  temperature,  concentration,  skin  friction  coefficient,  rate  of 
heat  transfer  and  rate  of  mass  transfer  at  the  plate  are  presented  in  graphs,  and  the  effects  of  various  physical  parameters 
like  Prandtl  number  Pr,  thermal  Crash  of  number  Gr,  mass  Crash  of  number  Cm,  Schmidt  number  Sc,  chemical 
reaction  parameter  Kr,  magnetic  field  parameter  M,  thermal  radiation  parameter  Nr,  Radiation  absorption  parameter  R 
and  heat  source  parameter  Qon  the  above  flow  quantities  are  analyzed  and  the  results  obtained  are  physically  interpreted. 
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INTRODUCTION 

The  phenomenon  of  mixed  convection  has  been  the  object  of  extensive  research.  The  importance  of  this 
phenomenon  is  increasing,  due  to  its  enhanced  concern  in  science  and  technology,  about  bounce  induced  motions 
in  the  atmosphere,  and  the  bodies  of  water  and  quasi  solid  bodies,  such  as  earth  are  related  to  them.  Convective 
flows  in  porous  media  have  applications  in  geothermal  energy  recovery,  oil  extraction,  thermal  energy  storage  and 
flow  through  filtering  devices.  The  hydro  magnetic  convection  with  heat  and  mass  transfer,  in  porous  medium  has 
been  studied.  It  is  due  to  its  importance  in  the  design  of  magneto  hydrodynamic  generators,  accelerators  in 
geophysics,  the  design  of  underground  water,  energy  storage  system,  astrophysics,  soil  sciences  and  so  on.  On 
account  of  their  importance,  these  flows  have  been  studied  by  several  authors,  among  them  are  Elbashbeshy  [1], 
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Yih  [2],  Hossain  and  Rees  [3],  Mahapatra  et  al.  [4],  Srinivas  and  Muthuraj  [5],  Raju  and  Varma  [6]. 

Combined  heat  and  mass  transfer  problems  with  chemical  reaction,  are  of  importance  in  many  processes  and 
have,  therefore,  received  a  considerable  amount  of  attention  in  recent  years.  In  the  processes  such  as  drying,  evaporation  at 
the  surface  of  a  water  body,  energy  transfer  in  a  wet  cooling  tower  and  the  flow  in  a  desert  cooler,  heat  and  mass  transfer 
occur  simultaneously.  Possible  applications  of  this  type  of  flow  can  be  found  in  many  industries.  For  example,  in  the  power 
industry,  among  the  methods  of  generating  electric  power  is  one  in  which,  electrical  energy  is  extracted  directly  from  a 
moving  conducting  fluid.  Chamkhaet  al  [7],  used  the  blotter  difference  method,  to  study  the  laminar  free  convection  flow 
of  air  past  a  semi-infinite  vertical  plate,  in  the  presence  of  chemical  species  concentration  and  thermal  radiation  effects. 
The  effect  of  chemical  reaction  to  free  convective  flow  and  mass  transfer,  of  a  viscous  incompressible  and  electrically 
conducting  fluid  over  a  stretching  sheet  was  investigated  by  Afify  [8],  in  the  presence  of  transverse  magnetic  field.  Ibrahim 
and  Abdullah  [9]  derived  the  analytical  solution  for  the  steady  MHD  mixed  convection,  laminar,  heat  and  mass  transfer 
flow,  over  an  isothermal  inclined  permeable  stretching  sheet,  immersed  in  a  uniform  porous  medium  in  the  presence  of 
chemical  reaction,  thermal  radiation,  Dufour  and  Soret  effects,  an  external  transverse  magnetic  field  and  internal  heating. 
In  the  study  of  Alam  et  al.  [10]  an  analysis  is  carried  out  to  investigate  the  effects  of  variable  chemical  reaction, 
thermophoresis,  temperature-dependent  viscosity  and  thermal  radiation  on  an  unsteady  MHD  free  convective  heat  and 
mass  transfer  flow  of  a  viscous,  incompressible,  electrically  conducting  fluid  past  an  impulsively  started  infinite  inclined 
porous  plate.  The  results  show  that,  higher  order  chemical  reaction  induces  the  concentration  of  the  particles  for  a 
destructive  reaction,  and  reduces  for  a  generative  reaction.  Rashad  et  al.  [11]  considered  MHD  free  convective  heat  and 
mass  transfer  of  a  chemically  reacting  fluid,  from  radiate  stretching  surface  embedded  in  a  saturated  porous  medium.  The 
Uri  et  al.  [12],  studied  unsteady  double  diffusive  MHD  boundary  layer  flow  of  chemically  reacting  fluid,  over  a  flat 
permeable  surface. 

The  study  of  heat  generation  or  absorption  effects,  in  moving  fluids  is  important  in  view  of  several  physical 
problems,  such  as  fluids  undergoing  exothermic  or  endothermic  chemical  reaction.  Due  to  the  fast  growth  of  electronic 
technology,  effective  cooling  of  electronic  equipment  has  become  warranted,  and  cooling  of  electronic  equipment  ranges 
from  individual  transistors  to  mainframe  computers  and  from  energy  suppliers  to  telephone  switchboards,  and  thermal 
diffusion  effect  has  been  utilized  for  isotope  separation  in  the  mixture,  between  gases  with  very  light  molecular  weight 
(hydrogen  and  helium)  and  medium  molecular  weight.  Hossain  et  al.  [13],  studied  the  problem  of  natural  convection  flow 
along  a  vertical  wavy  surface,  with  uniform  surface  temperature  in  the  presence  of  heat  generation/  absorption.  Hady  et  al. 
[14],  investigated  the  problem  of  free  convection  flow,  along  a  vertical  wavy  surface  embedded  in  electrically  conducting 
fluid  saturated  porous  media,  in  the  presence  of  internal  heat  generation  or  absorption  effect.  Chemical  reaction  and 
radiation  effects,  on  unsteady  MHD  free  convective  fluid  flow,  embedded  in  a  porous  medium  with  time  dependent  suction 
with  the  temperature  gradient  heat  source,  was  studied  by  Seshaiah  et  al.  [15].  Recently,  Balamurugan  et  al.  [16]  analyzed 
the  effects  of  chemical  reaction,  thermal  radiation  and  radiation  absorption  on  unsteady  double  diffusive  free  convection 
flow  of  Kuvshinski  fluid  past,  a  moving  porous  plate  with  heat  generation,  under  the  influence  of  a  uniform  transverse 
magnetic  field.  Das  et  al.  [17],  investigated  MHD  boundary  layer  slip  flow  and  heat  transfer  of  nano  fluid  past  a  vertical 
stretching  sheet,  with  non-uniform  heat  generation/  absorption. 

This  paper  aims  at  studying  the  effects  of  chemical  reaction  and  radiation  absorption,  on  unsteady  MHD  free 
convection  fluid  flow,  embedded  in  a  porous  medium  with  time  dependent  suction  and  temperature  gradient  heat  source,  in 
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two  different  cases  of  boundary  conditions. 


Case  (I):  u  —  u p ,  0  =  l  +  £ent,  C  =  l  +  £ent  at  y  =  0 


u  — ^  0,  0  — ^  0,  C  — ^  0  cis  y  — ^  °° 

Case  (II):  u  =  l  +  £ent,  0  =  l  +  £em,  C  =  l  +  £e'"  at  y  =  0 

u  — ^  0 ,  0  — ^  0 ,  C  — ^  0  as  y  — y 00 

FORMULATION  OF  THE  PROBLEM 

We  consider  a  problem  of  unsteady  two  dimensional,  laminar  boundary  layer  flow  of  viscous,  incompressible, 
electrically  conducting  fluid,  along  a  semi-infinite  vertical  plate  in  the  presence  of  thermal  and  concentration  boundary 

effects.  A  variable  time-dependent  suction  velocity  V  =  — V0(l  +  £Ae'U  )  is  considered  normal  to  the  How.  The  plate  is 

taken  along  x  -  axis  in  vertical  upward  direction,  against  the  gravitational  field,  y  -  axis  is  taken  normal  to  the  flow,  in  the 
direction  of  applied  transverse  magnetic  field.  Further,  due  to  semi-infinite  plane  surface  assumption,  the  flow  variables  are 
the  functions  of  y  and  t  only.  The  following  assumptions  are  made  in  the  governing  equations. 

•  The  viscous  dissipation  is  neglected. 

•  The  Joule  dissipation  is  neglected. 

•  The  induced  magnetic  field  is  assumed  to  be  negligible,  as  the  magnetic  Reynolds  number  of  the  flow  is  taken  to 
be  very  small. 

The  governing  equations  of  the  problem  in  dimensional  form  are 


(1) 


+  V 


*  du*  _  oBl 


u  +V 


(2) 


3/  p 


(3) 


— r  +  v  — r  =  v 
dt  dy 


dC  *  dC 


(4) 
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Figure  1:  Physical  Configuration  and  Coordinate  System 

The  corresponding  initial  and  boundary  conditions  are 

Case  (I):  u  =  u p  ,T  =  Tw  +  e(Tw  -rje" '  ,C  =  C  w  +  £(CW  -  Cx)e"  '  at  y- 0 

u*  — >  0,  T*  — >  7^  ,  C  — >  as  y  — >  °°  (5) 

Case  (II)  u=  \  +  SenV,  T  =  Tw  +  e{Tw  -  TJenY ,  C  =  Cw  +  e(Cw  -  CJenV  at  y  =  0 


u  — >0,  T*  — >7^,  C*  — >  C„  as  y— 


(6) 


By  using  Rosselant  approximation,  we  can  write  the  radiative  flux  £/,.  as: 


4<t*  3r*4 
3&J 


(7) 


It  is  assumed  that,  the  temperature  differences  within  the  flow  are  sufficiently  small  so  that  T  can  be  expanded 
in  a  Taylor  series,  about  the  free  stream  temperature  1 M  so  that,  after  rejecting  higher  order  terms: 


th  =4rY-3r4 

The  equation  of  energy  after  submission  of  equations  (7)  &  (8),  can  be  written  as 
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-TJ+R\C-CJ 

dt  dy  pcp  dy  3pcpk  dy*  pC p  dy 


(9) 


From  equation  (1),  v  —  —V0(l  +  FAe" '  )  (10) 

Where,  A  is  the  suction  parameter  and  £A«l.  Here,  the  minus  sign  indicates  that,  the  suction  is  towards  the 


plate. 


The  following  non-dimensional  parameters  are  introduced  in  the  equations  (2),  (4)  &  (9)  to  get  the  dimensionless 

form. 


U  y  U0 

u  =  —,  y  = - 

U0  V 


,jyi,  n=ii,  Pr=^A,  Sc=i,  m=ob"v 


u 


2’  "p  U'  k  D 


pu 


2  ’ 


Gr  = 


ygfiiK  -O 

ul 


k*U , 


,  k=—^,  e  = 


T  -T 


C-C 


1 6a  T: 


Gm  =  V^<-C~\  R  = 
Ul 


v“  Tv-T 

R*vK-cl) 


0  = — ; Nr  =  - 

C  -C  3k  k 


i  *2 

k,.  v 


pcATl-Tlwl 


’  Q- 


Qv 


P°  PUq 


(11) 


On  introducing  equation  (11)  into  equations  (2),  (4)  &  (9),  the  governing  equations  of  the  problem  in  non- 
dimensional  form  are 


dll 


nt  du  d2u 


dt  '  'dy  dy2 


- (1  +  eAen,)  —  =  —  +  GrO  +  Gm0 

dt 

do 


k  o 

M  +  — 

\  k  j 


u 


dt 

d0 

~dt 


(1  +  eAent)- 


do  _ 

fl  +  Nr 

dy 

l  Pr 

d0 

1  d20 

d20  „d0  , 

— v+Q 1-  R0 

dy2  dy 


-Kr2  0(14) 


dy  Sc  dy2 

The  corresponding  initial  and  boundary  conditions  are 

Case  (I):  u  —  u p ,  0  =  1  +  Fe"',  C-  l  +  £en‘  at  y-  0 

u  — y  0 ,  0  — y  0 ,  C  — y  0  cis  y  — y  °° 

Case  (II):  U  =  1  +  £e"  ,  0  =  1+Fe"' ,  C  =  1  +  fie'"  at  y  =  0 

u  — y  0 ,  0  — y  0 ,  C  — ^  0  as  y  — ^  °° 


(12) 


(13) 


(15) 


(16) 
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SOLUTION  OF  THE  PROBLEM 

The  equations  (12)  to  (14)  are  coupled  non-linear  partial  differential  equations  whose  solution,  in  closed  form  is 
difficult  to  obtain,  to  solve  these  non-linear  partial  differential  equations,  we  assumed  the  velocity,  temperature  and 
concentration  fields  as 

u  =  u0(y)  +  £en'ul(y)  +  0(£2) 

0  =  0o(y)  +  £en,0l(y)  +  0(£2)  (17) 

</>  =  </>0(y)  +  £e"t</> \(y)  +  0(£2) 

We  now  substitute  equation  (17)  in  equations  (12)  to  (14)  and  equating  harmonic  and  non-harmonic  terms. 


neglecting  higher  order  terms  in  £  ,  we  obtain: 

u0  +  u0  -MluQ  =  -Grd0  -  Gmt/>0  (18) 

Mj  +  Mj  -  ( n  +  M1)ul  =  ~Au0  -  Grd]  -  Gmtfiy  (19) 

0lt+h(\  +  Q)0o  -  -Rh</>0  (20) 

0t  +h(l  +  Q)9l  -  nh  9X  =  -M  0{)  -  Rh  (fly  (2 1 ) 

<p0  +  Sc <f)0  -  KrSctf>0  =  0  (22) 

(fly  +  Sc0x  -  Sc(Kr  +  n)(f)l  =  -ScA<f>0  (23) 


Where  M .  =  M  -I - ,  h  = - and  primes  indicate  differentiation  with  respect  to  y. 

k  1  +  Nr 

The  boundary  conditions  are: 

Case  (I):  U0=Up,  M,  =  0,  0Q  =  1,  <9,  =  1,  <pQ  =  1,  <j)x-  1  at  J  =  0 

i/q  — ^  0,  ux  — ^  0,  0q  — ^  0,  Oy  — ^  0,  (f)^  — ^  0,  (f)y  — ^  0  as  y  — ^  °°  (24) 

Case  (II):  M0  =  1,  Uy  —  1,  0O  =1,  <9,  =  1,  (f)(j  =  1,  (fly  =  1  y  =  0 

Uq  — ^  0,  Uy  — ^  0,  Oq  — ^  0,  ^  ^  0,  ^  0,  (fly  — ^  0  as  y  — ^  °°  (25) 

The  solutions  of  (18)  to  (23)  under  the  transformed  boundary  conditions  (24)  and  (25)  yield 

t=e~m'y  (26) 

(f)y  =Bye~nhy +B2e~m2y  (27) 
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0o=B3e~miy  +  (1-B3)e~m>y  (28) 

ex  =  B4e~m'y  +  B5e~m'y  +  B6e~m'y  +  57e“'"4V 
uQ  =  B%e~m'x  +  B,e~m'y  +  Bwe~nHy 

if  =  Bne~m'y  +  Bne~m'y  +  Bne~^y  +  Bue~m'y  +  Bl5e~"hy  +  Bx6e~m*y 

Skin  Friction 

The  skin-friction  coefficient  at  the  plate  is  given  by 


(29) 

(30) 

(31) 


(32) 


Nusselt  Number 


Rate  of  heat  transfer  in  terms  of  Nusselt  number,  at  the  plate  is  given  by 


Nu  = 


Sherwood  Number 


Rate  of  mass  transfer  in  terms  of  Sherwood  number  at  the  plate  is  given  by 


Sh  = 


v^Jy=0 


(33) 


(34) 


RESULTS  AND  DISCUSSIONS 

The  non-dimensional  boundary  layer  governing  equations  for  velocity  flow,  energy  and  concentrations  with 
corresponding  boundary  conditions  have  been  solved  by  perturbation  method.  In  the  present  study,  following  default 
parameter  values  are  adopted  for  computations  up=  0.5,  Sc  =  0.78,  t  =  1,  n  =  0.5,  Kr  =  0.5,  A  =  0.5,  e  =  0.02,  Nr  =  0.1,  Q  = 
0.1,  R  =  0. 1,  M  =  0.5,  k  =  1,  GM  =  5,  Gr  =  10.  All  graphs  therefore,  correspond  to  these  values  unless  specifically  indicated 
on  the  appropriate  graph. 

Case  (I):  The  fluid  velocity  variation  in  case  of  externally  cooled  plate  is  shown  in  figure  1  for  various  values  of 
Schmidt  number  Sc.  It  is  observed  that,  for  heavier  diffusing  foreign  species,  i.e.,  increasing  the  Schmidt  number  leads  to  a 
decrease  in  the  velocity.  Figure  2,  shows  the  effect  of  magnetic  parameter  M,  on  the  velocity  profiles.  The  velocity  profiles 
decrease  with  the  increase  in  M.  When  M  increases,  this  will  also  increase  the  Lorentz  force,  which  opposes  the  flow  and 
leads  to  enhanced  deceleration  of  the  velocity  profiles.  The  influence  of  the  radiation  absorption  parameter  R,  on  the 
velocity  and  temperature  profiles  is  shown  in  figures  3  and  4,  respectively.  It  is  obvious  that,  an  increase  in  the  radiation 
absorption  parameter  results  in  increasing  velocity  and  temperature,  within  the  boundary  layer. 

Case  (II):  Figure  5,  shows  that,  the  velocity  profiles  for  different  values  of  heat  source  parameter  Q.  The 
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numerical  results  show  that,  the  effect  of  increasing  values  of  heat  source  parameter  Q,  results  in  a  decreasing  velocity.  The 
influence  of  the  radiation  parameter  Nr  of  the  velocity  and  the  temperature  profiles  are  shown  in  figures  6  and  7, 
respectively. 


Figure  2:  Effect  of  Sc  on  Velocity  Field 


Figure  3:  Effect  of  M  on  Velocity  Field 
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Figure  4:  Effect  of  R  on  Velocity  Field 


Figure  5:  Effect  of  R  on  Temperature  Field 

The  radiation  parameter  Nr,  defines  the  relative  contribution  of  conduction  heat  transfer  to  thermal  radiation 
transfer.  It  is  obvious  that,  an  increase  in  the  radiation  parameter  results  in  increasing  velocity  and  decreasing  the 
temperature.  Figure  8,  shows  the  effect  of  chemical  reaction  parameter  Kr,  on  concentration  profiles.  It  is  observed  that,  an 
increase  in  chemical  reaction  parameter  Kr  decreases  in  concentration. 


www.tjprc.ore 


editor@tjprc.org 


Radiation  Absorption  Effects  on  Unsteady  MUD  Flow  of  Radiative  and  Chemically  Reacting 
Fluid  Past  a  Porous  Plate  with  Variable  Suction  and  Concentration 


171 


Figure  8:  Effect  of  Nr  on  Temperature  Field 


Figure  9:  Effect  of  Kr  on  Concentration  Field 

CONCLUSIONS 

The  problem  of  the  combined  influence  of  chemical  reaction  and  the  radiation  absorption  effect  of  a  magneto 
hydrodynamic  free  convective  flow,  embedded  in  a  porous  medium  with  time  dependent  suction  and  the  temperature 
gradient  heat  source,  has  been  analyzed  in  two  different  cases  of  boundary  conditions.  Graphical  results  for  various 
parametric  conditions  were  presented  and  discussed,  for  different  values.  The  main  conclusions  are  summarized  below. 

Case  (I): 

•  The  dimensionless  velocity  decreases  with  an  increase  in  Schmidt  number  search. 

•  As  the  magnetic  field  parameter  M  increases,  the  dimensionless  velocity  decreases. 
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•  With  the  increase  in  the  radiation  absorption,  the  velocity  temperature  increase. 

Case  (II): 

•  The  dimensionless  velocity  decreases  with  an  increase  in  heat  source  parameter  Q. 

•  Velocity  increases  with  an  increase  in  radiation  parameter  Nr. 

•  Temperature  decreases  with  an  increase  in  radiation  parameter  Nr. 

•  Concentration  reduces  with  an  increase  in  chemical  reaction  parameter  Kr. 
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